УДК 681.322 


Д.А. Вятченин, В.В. Старовойтов 
Объединенный институт проблем информатики НАН Беларуси, г. Минск 
уайсвепт(@лта!.га 


Идентификация объектов посредством 
системы нечеткого вывода 


В статье предлагается методология идентификации объектов на аэрокосмических изображениях в 
режиме времени, близком к реальному, посредством механизма нечеткого вывода на примере системы 
нечеткого вывода типа Мамдани. Кратко излагается метод быстрого прототипирования системы 
нечеткого вывода на основе результатов обработки данных об объектах обучающей выборки эвристическим 
алгоритмом возможностной кластеризации. Приводятся результаты вычислительных экспериментов 
и сформулированы предварительные выводы. 


Введение 


В задачах обработки и анализа аэрокосмических изображений особое место 
занимают методы идентификации обнаруженных на снимках объектов, то есть методы 
отнесения того или иного объекта к одному из нескольких заданных классов. Следует 
указать, что в ряде задач, таких, к примеру, как поддержание в актуальном состоя- 
нии кадастровых карт, мониторинге земной поверхности в целях обнаружения чрез- 
вычайных ситуаций, а также при разработке систем специального назначения, обнару- 
жение и идентификация объектов заданных классов должна производиться в реальном 
режиме времени. Для решения задачи идентификации объектов традиционно исполь- 
зуются, с одной стороны, нейросетевые модели, а с другой — продукционные модели. 

В условиях неточной, нечеткой и неопределенной информации представляется 
целесообразным обращение к нечетким моделям, что позволяет получить, с одной сторо- 
ны, точные, а с другой — содержательно осмысленные результаты идентификации объектов. 

Целью данной работы является разработка общей схемы применения системы 
нечеткого вывода типа Мамдани [1], генерируемой на основе результатов обработки 
данных об объектах обучающей выборки О-АЕС(с)-алгоритмом возможностной 
кластеризации [2], предусматривающей возможность дообучения соответствующей 
системы в процессе распознавания объектов на аэрокосмических снимках. 


Методы построения базы правил системы нечеткого 
вывода, основанные на нечеткой кластеризации 


Системы нечеткого вывода представляют собой алгоритм получения нечетких 
заключений на основе нечетких предпосылок [3]. Главным элементом систем нечет- 
кого вывода является база нечетких продукционных правил, в наиболее общем случае 
имеющих вид 


[: ЕСЛИ есть В!И...И "есть В” ТО у, есть СИ... И у, есть С", (1) 


где /е1....,с} — номер правила, & е №', 1Е1,....т\ — входные переменные и №’ — 
область определения соответствующей переменной, у ЕТ, /[Е{,...,с} — нечеткие 
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выходные переменные, причем у, — область определения соответствующего 


1 
заклюючения, а В,, С, — нечеткие множества с функциями принадлежности у ы (Е) 
1 


и Ус (у,), определенные на соответствующих универсумах. В правиле вида (1) т 


предпосылкам соответствует с заключений, в силу чего структура такого правила в 
специальной литературе именуется МПМО-структурой [3]. 

База нечетких правил, являющаяся основным элементом системы нечеткого 
вывода, может формироваться, с одной стороны, экспертным путем, а с другой - на 
основе обработки данных обучающей выборки, для чего чаще всего используются 
оптимизационные методы нечеткой или возможностной кластеризации [4]. При вто- 
ром подходе исследуемая совокупность Х = {х,,...,х,} объектов обучающей выборки 


обрабатывается каким-либо методом нечеткой или возможностной кластеризации с 
последующим проецированием значений принадлежности д,, или значений типи- 


1 
чности и, того или иного нечеткого кластера А’, [=1,....с на координатные оси 
признакового пространства /”(Х). Полученные таким образом нечеткие множества 


с дискретными функциями принадлежности У г. (©) аппроксимируются параметрически- 
1 
ми непрерывными функциями принадлежности У, (2) — как правило, треугольной 
Г 


или трапецевидной формы — для чего целесообразно воспользоваться алгоритмом, 
предложенным в [5] М. Сугэно и Т. Ясукавой. 

Так как в оптимизационных методах нечеткой и возможностной кластеризации 
первоначальное разбиение формируется случайным образом, то для получения при- 
емлемого результата классификации в виде нечеткого с -разбиения или возможностного 
разбиения, необходимым является проведение серии вычислительных экспериментов, 
что значительно затрудняет, или делает вообще невозможным генерирование базы 
нечетких правил в реальном режиме времени, или, по крайней мере, в режиме, близ- 
ком к реальному. Указанного недостатка лишен предложенный в [6] метод проекти- 
рования систем нечеткого вывода, основанный на обработке данных об объектах 
обучающей выборки эвристическим О-АЕС(с)-алгоритмом возможностной класте- 
ризации [2], результатом работы которого является так называемое распределение 
В*(Х)= А ..А„} объектов совокупности Х’по заданному числу с частично или 
полностью разделенных нечетких с -кластеров, которое представляет собой частный 
случай возможностного разбиения [7], а также соответствующее значение порога 
сходства а. Основные понятия и схема О-АЕС(с)-алгоритма изложены в [2], [8], [9]. 

Предложенный в [6] метод извлечения нечетких правил, основанный на 
обработке данных обучающей выборки, представленных матрицей «объект-признак» 


Хи =[8], 1=1...п, #=1,...т, где п - число объектов классифицируемого 
множества АХ’ = {х,...,х,}, а т — размерность признакового пространства, с помощью 
2-АЕС(с)-алгоритма, заключается в следующем. Нечеткие множества В+, [Е {1,...,с}, 
1Е{|,...,т} предпосылок продукционных правил вида (1) с, в общем случае, 
трапецевидными функциями принадлежности у,„(&) могут быть представлены в 
1 
Я 1 ЕЕ 5 1 
параметрической форме В, = (а„,т,„, ти, а), где т\, — нижнее модальное значение 
нечеткого интервала В’, т, - его верхнее модальное значение, а величины 


ГА 1 — — > 
т„-ар) и (а»)-т,) представляют собой значения левого и правого коэффи- 


«Штучний 1нтелект» 32010 313 


—4в_ Вятченин Д.А., Старовойтов В.В. 


циентов нечеткости соответственно [10]. Необходимо отметить, что в случае 

совпадения нижнего и верхнего модальных значений, т, = =т,, нечеткий 

интервал В, предпосылки продукционного правила, изображенный на рис. 1 а), 
> 1 Е 

будет представлять собой треугольное нечеткое число В, =(а,,ти,@и) [10]. При 


построении условных частей нечетких правил задача заключается в определении со- 
ответствующих параметров нечетких множеств В,, [Е\|...,с}, [Е \1,...,т}. 


те, 


1=1,....И, #=1,....т, заключающейся в нормировке значений %, к примеру, с 


После предварительной обработки матрицы исходных данных Х 


пхт 


помощью нормализации 


х=, (2) 


так что каждая строка х, матрицы Х’„„„ =[х!], [=1,....И, #=1,....т может интер- 
претироваться как нечеткое множество на универсуме признаков с соответствующей 
функцией принадлежности х, = ты (х’) = [01], 1=1.....п, #=1,....т, и применением 
к полученному множеству и нечетких множеств какого-либо расстояния между 


нечеткими множествами [11] — например, квадрата нормализованного евклидова 
расстояния 


т 


обеих = и), 6} (3) 


15 
получается матрица нечеткого отношения несходства /=[м,(х,‚,х,)], 5 Л=1....,П, 
применяя к которой, в свою очередь, операцию дополнения 
Из (хьх,)=1-и,(х,ьх,), МЕ И (4) 
строится матрица нечеткой толерантности Т’=[4,(х,,х,)], Ь.] =1....п, на универсуме 
над 
Для построения базы правил системы нечеткого вывода типа Мамдани [1] 


данные об объектах исследуемой совокупности, представленные в форме матрицы 
нечеткой толерантности Т=[м,(х,,х,)], Ь/=1....п обрабатываются О-АЕС(с)- 
алгоритмом для заданного числа классов с, и для каждого нечеткого а -кластера 


1 * 1 с 
А’ „› [Е {1,...,с} полученного распределения К”(Х)={А„,,..., А„)} выделяется его 


ядро К (А) и носитель А’. Кроме того, как указывалось выше, дополнительным 
результатом классификации является значение порога сходства а е (01|, соответ- 
ствующего полученному результату классификации А”(Х). Интервал [из › Жиних | 


значений каждого признака ®, (Е 11,....т} вычисляется для носителя ри каждого 


1 
(а) › 


может быть получено по формуле 
& = ши, У!Е1,....т\, У! Е 1... с}, (5) 


а 


нечеткого @-кластера А 1Е\,...,с}. В частности, значение Е 1Е\|,...,т} 


Л) шм 
1Е\,...,т} — по формуле 
& = тах, У!Е\,...,т}, УГ Е 1,...,с}. (6) 
ХЕ А, 


(Г) тах 
© 


а значение * 


(Г) тах › 
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В свою очередь, параметр а вычисляется путем решения линейных уравнений 


при условиях 


Тв (бит) = 1 - @), Ив (@)=0, (7) 
а параметр а) — при условиях 
в (раю) = (1-а), в (Я) =0. (8) 


Значение Е вычисляется для всех типичных точек г, ЕК (Аш) нечеткого &- 


кластера и ‚ [Е {|,...,с} в соответствии с выражением 


Е ‚0 #, Че>\....,й}, (9) 
ТЕК (Аа)) 
а значение *, — по формуле 
„= шах ®, Уее{....,/}. (10) 


пеЕК(Аа)) 
1 > 
Таким образом, параметр т, может быть вычислен путем решения линейных 


уравнений, исходя из условий 


Ув (и )= У: (ти) =1, (11) 
а параметр т) вычисляется аналогичным образом при условиях 
Уи (Во) =Ун (ти) =1. (12) 


Таким образом, число входных переменных & генерируемых правил вида (1) 


т 


соответствует размерности признакового пространства /”(Х), а число выходных 


переменных у,, которым в соответствие могут быть поставлены метки классов, как 
и число продукционных правил, эквивалентно числу нечетких @ -кластеров Я 
1=Е11,...,с} полученного распределения К“(Х). 


Функции принадлежности заключений нечетких правил целесообразно строить 
следующим образом: нечеткие множества С/, [Е1....,с} заключений продукцион- 


ных правил вида (1) представляются в параметрической форме С! = (а, и „№ ), где 
@ — порог сходства, при котором строится решение А“”(Х) задачи классификации, 


также являющийся, как указывалось выше, результатом работы О-АЕС(с)-алгоритма, 
ш, - наименьшее значение принадлежностей объектов х, е А' соответствующему 


нечеткому & -кластеру и ЕА”(Х), и, = и и, а и, — наибольшее значение принад- 
Е* ХЕ Ад 


лежности объектов хе А’ нечеткому @-кластеру А) =А"(Х), 1Е\....,с}, 


и, = тах д, ‚ причем, если в результате обработки данных О-АЕС(с)-алгоритмом полу- 


7 а 


ченные нечеткие & -кластеры Яо ЕА”(Х), [Е1,...,с} будут являться сильными не- 

четкими & -кластерами [12], то, очевидно, будет иметь место д, =1, так что трапе- 

цевидная функция принадлежности 7, (у,) соответствующего нечеткого множества 
1 

С', [Е1....,с} примет вид, изображенный на рис. 1 6). В случае же, когда носитель 


1 * 
А, некоторого нечеткого @ -кластера Я построенного распределения А”(Х) пред- 
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ставляет собой пустое множество, то выходной переменной у,, [е {1,...,с} будет соот- 


1 
ветствовать пустое нечеткое множество С, , так что У„„(у,)=0. 
1 


7, (8) 7е (у) 


гл 


у 1 — м1 ее 
бут 71. 


1 я 
о Житак & ты Щ=1 


Рисунок 1 — Функции принадлежности нечеткого множества: 
а) предпосылки; 6) заключения 


Правила построенной указанным способом базы нечетких продукций будут 
иметь вид (1) как в случае, когда нечеткие а -кластеры и 1Е\|,....с} являются 
частично разделенными, так и в случае полностью разделенных нечетких а -клас- 


теров, причем в случае частично разделенных нечетких а -кластеров, каждой из с 
условных частей правил (1) будут соответствовать несколько непустых нечетких 


множеств С’, [Е 1|,...,с}, соответствующих выходным переменным у’, [Е {1....,С}. 


Формирование обучающих данных для построения 
системы нечеткого вывода 


Для построения системы нечеткого вывода необходимо подготовить обуча- 
ющую информацию, характеризующую объекты соответствующих классов. Так как 
распознаванию подлежат присутствующие на снимках объекты, то процесс подго- 
товки данных для последующей обработки О-АЕС(с)-алгоритмом с целью постро- 
ения системы нечеткого вывода может быть в виде последовательности этапов, 
представленных в табл. 1. 


Таблица 1 — Формирование данных для построения системы нечеткого вывода 


Номер и содержание этапа Методы решения соответствующей задачи 

1 | Предварительная обработка | Контрастирование, фильтрация и т.п. 
изображения 

2 | Локализация области Интерактивное указание оператором; 
интереса, содержащего автоматическое определение 
объект определенного класса 

3 | Обработка фрагмента Сегментация изображения 
изображения 

4 | Выделение объекта на Устранение шума; удаление областей, не 
изображении относящимся к объектам интереса 

5 | Вычисление характеристик В зависимости от задачи, конкретизируются 
обнаруженного объекта геометрические количественные 

характеристики: длина, периметр и т.д. 
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При вычислении характеристик объектов следует учитывать, что объекты одного 
класса могут ориентироваться на изображении различным образом, так что характе- 
ристики объектов ДОЛЖНЫ выбираться с учетом данного обстоятельства. К примеру, 
характеристиками объекта могут служить: длина в пикселях проекции объекта на ось 
Х, длина проекции на ось У, а также общее число пикселей, составляющих объект 


интереса, как это изображено на рис. 2. 


. 


ь 


ЯР 


т 


Нах к 


Рисунок 2 — Признаки объекта 


Упрощенный пример процесса обнару- 
жения объекта заданного класса на исходном 
изображении приведен на рис. 3 -— так, на 
рис. 3 а) представлено исходное изображение 
размерностью 1277х928 пикселей, на рис. 
3 6) — интерактивно выделенный фрагмент 
исходного изображения размерностью 197х183 
пикселя, предположительно содержащий объ- 
ект класса «автомобили», на рис. 3 в) изоб- 
ражен результат сегментации выделенного 
фрагмента, представляющий собой бинарное 
изображение, содержащее как связную об- 
ласть, возможно, являющуюся объектом задан- 
ного класса, так и другие связные области, 


и, наконец, на рис. 3 г) представлен результат после удаления шума и исключения 
других связных объектов — бинарное изображение, определяющее расположение 


объекта на снимке. 


6) 


г) 


Рисунок 3 — Этапы обнаружения объекта заданного класса на изображении: 
а) исходное изображение; 6) фрагмент изображения; в) бинарное изображение — 
результат сегментации фрагмента; г) бинарное изображение объекта на фоне 
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Следует отметить, что для сегментации изображений могут использоваться самые 
разнообразные подходы, такие, как методы кластеризации, метод управляемого 
водораздела, а также пороговая сегментация [13]; при этом необходимо указать, что 
выбор того или иного метода сегментации диктуется условиями конкретной задачи и 
зависит от особенностей обрабатываемого изображения. 


Общая схема применения системы нечеткого вывода 
для идентификации объектов на изображениях 


При идентификации с помощью системы нечеткого вывода вновь поступающих 
объектов может оказаться, что объект не будет отнесен ни к одному из заданных 
классов. В такой ситуации система нечеткого вывода должна быть сгенерирована 
повторно с учетом данных о нераспознанном объекте. Общая схема процесса генери- 
рования системы нечеткого вывода с учетом данных о поступающих на класссифи- 
кацию объектов изображена на рис. 4. 


Данные об объекте, 
поступающем на 
распознавание 


Данные об объектах 


обучающей выборки 


Определение числа Объект 
классов в обучающей идентифицирован? 
выборке 


Еластеризация данных с помощью 


Р 
Т-АРС/с)-алгоритма и построение Е 


распознавания 


СИСТЕМЫ НЕЧЕТКОГО ВЫВОДА 


Рисунок 4 — Общая схема дообучения системы нечеткого вывода 


На втором шаге представленной выше схемы для определения числа классов, 
на которое «расслаивается» обучающая выборка, можно использовать любую быст- 
рую кластер-процедуру, не требующую априорного задания числа классов с -— к 
примеру, О-АЕС-ТС-алгоритм возможностной кластеризации [14] или субтрак- 
тивный 5СМ-алгоритм [15]. Следует, однако, учитывать, что при использовании О- 
АЕС-ТС-алгоритма вид нормировки исходных данных и вид функции расстояния 
должен совпадать с нормировкой и расстоянием, используемыми при кластеризации 
данных О-АЕС(с)-алгоритмом в процессе построения системы нечеткого вывода. 


Иллюстративный пример 


Представляется целесообразным проиллюстрировать сущность предложенной 
методологии на простом примере. Для проведения вычислительных экспериментов 
были выбраны 50 фрагментов изображений, на которых представлены 30 объектов 
класса «автомобили» и 30 объектов, не относящихся к указанному классу, такие, как 
элементы дорожной инфраструктуры, бассейны, другие сооружения и т.п. Сегмен- 
тация изображений была выполнена с помощью ЗСЕСМ-алгоритма нечеткой класте- 
ризации [16] с последующим выделением а -ядер кластеров полученного нечеткого 
с-разбиения [17], причем в качестве признакового пространства была выбрана 
яркость в каждой из компонент цветового пространства КОВ. 
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В качестве признаков распознаваемых объектов были выбраны рассмотренные 
выше: АХ - длина в пикселях проекции объекта на ось Х, АУ — длина проекции на 
ось У, а также МитРих — общее число пикселей, составляющих объект. Исследуемая 
совокупность была разделена на 2 части — обучающую выборку, включающую изо- 
бражения 20 автомобилей и 20 других объектов, и тестовую, включающую изобра- 
жения 10 автомобилей и 10 других объектов. Проекция выбранных для проведения 
экспериментов данных на плоскость первых двух признаков представлена на рис. 5, 
где объекты обучающей выборки обозначены символом о, а объекты тестовой — 
символом м. 
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Рисунок 5 — Исходные данные для вычислительных экспериментов 


В первом эксперименте с помощью О-АЕС-ТС-алгоритма [14] возможностной 
кластеризации с использованием нормализации (2), квадрата нормализованного ев- 
клидова расстояния (3) и дополнения (4) были обработаны данные об объектах обучаю- 
щей выборки, в результате чего число классов, на которые «расслаивается» указан- 
ная совокупность, оказалось равным 3. При разбиении данных О-АЕС(с)-алгоритмом 
на 3 класса, также с использованием формул (2), (3) и (4), оказалось, что 4 объекта 
оказались классифицированными некорректно — объекты, не являющиеся автомоби- 
лями, были классифицированы как автомобили. Распределение К”(Х) по 3 пол- 


ностью разделенным нечетким кластерам было получено при значении порога сход- 
ства а = 0,9347. Исходя из результатов кластеризации, была построена система 
нечеткого вывода, включающая 3 правила. На рис. 6 приведен пример работы по- 
строенной системы при идентификации объекта, изображенного на рис. 3 г). 
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Рисунок 6 — Пример классификации объекта исследуемой совокупности с помощью 
построенной системы нечеткого вывода типа Мамдани 
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В процессе распознавания всех 60 объектов некорректно идентифицирован- 
ными оказались 4 объекта, которые также были некорректно классифицированными 
при обработке обучающей выборки О-АЕС(с)-алгоритмом, и 7 объектов, являю- 
щихся элементами тестовой выборки, оказались неидентифицированными. Таким 
образом, из 20 объектов — элементов тестовой выборки, 13 объектов оказались клас- 
сифицированы корректно. Следует также отметить, что среди неидентифицирован- 
ных объектов не оказалось объектов класса «автомобили». 

Во втором эксперименте данные об объектах тестовой выборки были добав- 
лены к данным об объектах обучающей выборки, и с помощью О-АЕС-ТС-алгоритма 
было определено, что вся совокупность также «расслаивается» на 3 класса. При 
кластеризации данных о всей совокупности с помощью О-АЕС(с)-алгоритма, рас- 


пределение А“"(Х) по 3 полностью разделенным нечетким кластерам было получено 
при а = 0,9402 , причем некорректно классифицированными оказались 4 объекта, как 


и в первом эксперименте. На основании результатов кластеризации была сгенериро- 
вана система нечеткого вывода, и в процессе распознавания некорректно иденти- 
фицированными оказались 4 объекта, которые были некорректно классифицирован- 
ными О-АЕС(с)-алгоритмом, а неидентифицированных объектов не оказалось. Таким 
образом, на простом примере продемонстрирована эффективность предложенной 
методики дообучения системы нечеткого вывода. 


Заключение 


Предложенная методология построения системы нечеткого вывода позволяет 
генерировать систему нечеткого вывода в режиме времени, близком к реальному; 
кроме того, предложенный подход позволяет дообучать систему на основе данных о 
неидентифицированных объектах. Помимо точности и релевантности результатов 
классификации объектов, полученных с использованием построенной с помощью П- 
АЕС(с)-алгоритма возможностной кластеризации системы нечеткого вывода, глав- 
ным достоинством предложенного подхода к проектированию системы является 
возможность обработки данных в полностью автоматическом режиме. 
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Д.А. Вятчент, В.В. Старовойтов 

Гдентиф!кащя об’ектв за допомогою системи неч#ткого висновку 

У статт! пропонуеться методоломя 1дентифлкаци об’ект!в аерокосмачних зображеннях у режим! часу, 
близькому до реального, за допомогою механ1зму нечткого висновку на приклад! системи нечткого 
висновку типу Мамдан!. Стисло викладаеться метод швидкого прототипування системи нечиткого 
висновку на основ! результатв обробки даних об об’ектах навчально! виробки евристичним 
аморитмом можлив!стно! кластеризаци. Наводяться результати обчислювальних експериментв та 
сформульован! попередн! висновки. 


р.А. Иайисйетт, И.И. багоуойоу 

ОЦес5 14епиЙсайоп ТЬгопов Еи7ту Пегепсе Зузет 

ТБе рарег ргорозез а те#о4о|огу ю 1епИБу Фе обес оп аегозрасе ппазез ш-@те, с1озе 1ю Ще геа1, 
ФгоцэЬ бе тесБап15т оЁ Ра7ху шЕегепсе Юг Фе ехатр!е оГ Матдапг’ $ фуре Ёа7ху шЕегепсе зузет. Тре 
тефо4 оЁ гар! ргоюбурше оЁР а №7ту шЁЕегепсе зубет Базе оп йе тез оЁ Фе Наше 4аёа зе 
ргосеззте Бу Вейлизис а|еогИбт роз$ЬШзИс сшазегте Ш Чезстфед ЪФиеНу. Кези оЁ питепйса1 
ехрегитеп! аге ргезеще4 ап4 ргешитагу сопс[аз1оп$ аге Фогии!е4. 
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